Cooling of a mirror in cavity optomechanics with a chirped pulse 
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We investigate the response of a harmonically confined mirror to an optical pulse in cavity op- 
tomechanics. We show that when the pulsed coupling strength takes the form of a chirped pulse, 
thermal fluctuations of the mirror can be significantly transferred to the cavity field. In addition, 
the frequency modulation of the pulse could enable a better cooling performance by suppressing the 
sensitivity of the dependence of detuning and pulse areas. Using numerical investigations, we find 
that the pulsed-cooling is mainly limited by the cavity-field decay rate. 
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I. INTRODUCTION 

In cavity optomechanical systems [H-[l] , the cooling of a 
mechanical resonator is important to study the mechan- 
ical effects of light in the quantum domain [El-l22j. In 
addition, the reduction of thermal noise generally is a re- 
quirement to implement various applications in quantum 
information relying on optomechanical couplings [23U26I ] . 
With the rece nt progr ess of cooling in cavity optomechan- 
ical systems [27H36| . it is becoming possible to access 
quantum ground states. The cooling techniques, such as 
the resolved sideband cooling [27H291 [3J-|36|], generally 
make use of a continuous light field that drives a me- 
chanical resonator into a steady state of lower temper- 
ature. Recently, several authors have begun to explore 
the possibilities of cooling and manipulating the states of 
mirrors with optical pulses [3^ - l4oT | . 

Motivated by the fact that chirped pulses can lead to 
an efficient population transfer in two-level systems (4ll — 
[46j, we examine how a chirped-pulse interaction can 
transfer thermal fluctuations from the mirror to the cav- 
ity field. As we shall discuss below, the correspondence 
between a two-level system and optomechanical systems 
[defined in Eq. ([1])] can be established via the Heisenberg 
equations of motion of the linearized system under a ro- 
tating wave approximation (RWA). Therefore, analytical 
solutions known in two-level systems driven by chirped 
pulses can be applied here [46| . Owing to the frequency 
modulation in chirped pulses, population transfer can be 
made without the need for a high-precision control of 
detuning and pulse areas. 

In this paper, we will first provide a formulation of 
the linearized quantum system driven by a general pulse. 
Then we will indicate how the mirror-field coupling can 
be shaped into a chirped form by using a proper time- 
dependence of an external driving field. Specifically, we 
will study a class of chirped pulses proposed by Allen and 
Eberly for two-level systems [46[ . Such a class of chirped 
pulses has analytic solutions for nondissipative systems, 
and depending on the chirped parameters, they describe 
both adiabatic and nonadiabatic transitions. Using nu- 
merical calculations, we include counter-rotating terms 
and dissipation effects, and we demonstrate that cooling 



can be achieved with the chirped-pulse coupling. In the 
case where the mechanical damping rate of the mirror 
is sufficiently small, the cooling performance is mainly 
limited by the cavity-field decay rate. 



II. MODEL 

Our model consists of a Fabry-Perot cavity formed by 
a fixed end mirror and a moving end mirror connected 
with a spring (Fig. [l}. We consider a single cavity field 
mode with a resonance frequency lj c and creation (anni- 
hilation) operator (a). The moving mirror is treated 
as a quantum harmonic oscillator with a frequency u m 
and creation (annihilation) operator (b). Assuming 
the cavity is driven by an external field with a carrier 
frequency wi and a time- varying amplitude f2(i), the 
Hamiltonian of the system (in a rotating frame with the 
frequency lot) is given by, 

Hg — HA c d'a + fiLJ m b'b — hga'aQr+b) 

+Ml(t)(J + Ml*(t)a, (1) 

where A c = w c -W£ is the detuning and g is the radiation- 
pressure coupling strength. 

To treat the damping and noise in our model, we 
consider the system linearly coupled to oscillator baths. 
Under the Markovian approximation and neglecting 
counter-rotating terms in system-bath coupling, the 
quantum Langevin equations for the operators a and b 
are given by, 

a = —iA c a + iga(b* + b) — ifl(t) + dj n , (2a) 



b = —iu m b + iga) a 



-b + h 



(2b) 



where 7 C ("f m ) is the cavity field (mirror motion) de- 
cay rate, aj„ is the vacuum radiation noise operator 
for the cavity and 6j„ is the mechanical noise operator 
for the mirror. Both a,i n and fcj n have zero mean val- 
ues and they are characterized by the correlation func- 
tions {a ln (t)a\ n (t')) = lc 5(t-t'), (a\ n {t)a m {t')) = 0, 
(bin{t)bt n {t')) = lm (n m + l)8(t-t'), and {b\ n {t)b ln {t')) = 
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FIG. 1: (Color online) Schematic diagram of the cavity op- 
tomechanical system. A Fabry-Perot cavity, formed by a fixed 
end mirror and a harmonically bound end mirror, is driven 
by a pulse. 

Jmn m 8(t-t'), where n m = [exp(fiw m /fc B T' m )-l] _1 is the 
average thermal excitation number of the mirror at tem- 
perature T m . In this paper we will investigate the regime 
with uj rn 3> 7m- This is a regime where the Markovian 
approximation for the mirror noise can be justified [l8j . 



with A(t) = A c - 2gRe[(b(t)}]. The (a(t)) and (b(t)) are 
governed by: 

(a) = -iA(t)(a)-ia(t)-^(a), (4a) 

(b) = -iu; m (b)+ig\(a)\ 2 -^(b). (4b) 

For the linear approximation made above, we have ne- 
glected nonlinear terms igSa(5b + 5b^) in Eq. (|3aj) and 
igSaUa in Eg. (f3b| 

Equation ([3]) corresponds to a linear coupling described 
by the Hamiltonian iJ 7 = -%[(a t (*))5a+(a(t))<Ja t ](<S& 1 '+ 
5b). Here (a(t)) modulates the mirror- field coupling and 
its time dependence can be controlled by the driving am- 
plitude f2(i). We point out that any desirable (a(t)) as a 
function of time can in principle be achieved by a corre- 
sponding ft(t), according to Eq. (4a). For convenience, 
we let 



III. LINEARIZED SYSTEM AND FORMAL 
SOLUTION 

By writing o = (o) + 5o (o = a, b) and assuming the 
fluctuations are small (|(o)| 2 ^ (So^So)) during the pulse 
interaction, we may linearize Eq. ^) and obtain the equa- 
tions of motion, 

5a = —iA(t)5a + ig{a)(5b^ + 5b) — -^5a + a, n , (3a) 
5b = —ioj m 5b + ig[{a))5a + (a)5a < ] — ~^fib + bi n , (3b) 



g(a1(t)) = x (t)e l ^e- 2l9 ti Rc ^ T ^ dT , (5) 



where x(t) and 4>(t) are real functions, and the phase an- 
gle — 2g J* Re[(6(r))]<ir is introduced in order to compen- 
sate for the phase shift induced by the dynamical cavity 
frequency shift in A(t). 

By defining operators 5A(t) = 5ae^ {t)+f « A(r H 
and SB (t) — 8be lulmt , Eq. ([3]) can be concisely writ- 
ten as v(t) = M(t)v(i) + N(t), where v(t) = 
[5A(t),5B(t),5A^(t),5B i <(t)} T , and 



■ 2a +*^(t) 



M(t) 



ix(i)e l(A ^" m)t 
ix(i)e" 4(Ac ^ m)t ix(i)e l(Ac+Wm) * 
-i X (t)e- l ^+^ t _^-i0(f) 



-ix(t)e 



-i(A c +uj m )t 







-ix(i)e 



(A c -ui m )t 





-zxWe" l(Ac " Wm)t 



(6) 



and N(t) = [A m (t),B m (t),Al(t),Bl(t)] T with 

A in (i) = a^e^*^ A «<H and B^i) = 6 m e^*. The 
solution of v(t) is 

v(i) = G(t)v(0) + G(t) f G-^rjNCrJdr, (7) 
Jo 

where G(t) is governed by 

G(i) = M(t)G(f), (8) 
with G(0) = I being the identity matrix. 



The state of the system can be conveniently de- 
scribed by a covariance matrix R(f) whose elements 
are: R w (t) = (vj(t)vj/(t)) (1,1' = 1,2,3,4). Therefore 
R31M = (Sa^Sa) and R42(i) = (Sb^Sb) are mean dis- 
placed particle numbers measuring the fluctuations. By 
Eq. (7), R(t) reads, 

R(i) = G(i)R(0)G T (i) + G(t)Z(i)G T (i), (9) 

with 

z(t)= / / G-V)C(Ty)[G- 1 (<r / )] T *-dT'. (10) 

Jo Jo 
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The matrix R(0) is determined by the initial condition of 
the system, and C(t, t') is the two-time correlation func- 
tion of noise operators, which is defined by the elements 
C m ,(t,t') = <N z (t)N,,(t')> (/,/' = 1,2,3,4). Assuming 
that initially the cavity is in vacuum and the mirror is 
in a thermal equilibrium at the same temperature T m as 
its bath, i.e., p(Q) = |0) c (0| c <g> Pth {T m ) with Pth {T m ) = 
exjp[-hkj m b*b/kBT m ]/Tr(exp[-hkj m b*b/k B T m ]), then 
the matrix R(0) has three nonzero elements: Ri3(0) = 1, 
R24(0) = n m + 1, and R42(0) = n m . In addition, the 
Markovian baths imply C(r, r') = C5(t — r'), where 
C is a constant matrix with three nonzero elements: 
Ci 3 = 7c C 24 = J m (n m + 1), and C 42 = J m n m . 

IV. COOLING OF THE MIRROR 




CO t Q) t 

m m 

FIG. 2: (Color online) Plot of (a) the chirped pulse given 
in Eq. ()lip and (b) the original driving amplitude Q(t) vs 
the scaled time u) m t. The parameters are: ot/u m = 0.14, 
/3/u m = 0.04, uj m t = 40, and X o = |a/q 2 + P 2 - 



We will employ the expectation value of displaced 
phonon number {5b* 8b) as an indicator of cooling. This 
means that the idea of cooling in our scheme should be 
understood as a process of reducing excitations with re- 
spect to the mean amplitude (b) of the mirror. 



A. Chirped-pulse coupling 

We now ask what il(t) is suitable for cooling the mir- 
ror. Guided by the fact that a chirped-pulse driving 
can efficiently realize population transfer in two-level sys- 
tems [4H - I46I ] , we consider the coupling strength in Eq. ([5]) 
taking the chirped form (4o| . 

X(t) = Xosech[a(i - t )], (Ha) 
0(f) = <j)(t) = /3tanh[a(t -t )]. (lib) 

Here to determines the time of the pulse peak entering 
the cavity, a -1 measures the pulse duration, (3 controls 
the magnitude of the frequency modulation, and xo is 
the strength of the pulse coupling. 

The required driving amplitude Q,(t) for generating the 
above x(t) and <f>{t) can be found using Eqs. flU, (0, and 
(gUl as 



fi(i) = i{a{t)) - A(t) -C^ (o(t)>, 



(12) 



with 



(a(t)) 



Xo 
fJ 



sech[a(t - to ) ]e -i[m-2 9 ft**KHr))]dT} 



(13a) 



(b(t))=i^ / sech 2 [a{r -t )]e-^ + ^) {t - T) d T , 
9 Jo 

(13b) 



log 



cosh[a(t — t )] 



cosh(a£o) 



(13c) 



An example illustrating the chirped pulse coupling and 
the corresponding f2(t) is given in Fig. [5] 



B. Ideal case 

To find the optimal relation for xo, 01, and (3 in Eq. (jlip 
such that efficient cooling can be achieved, we first in- 
vestigate the nondissipative case (j c = and j m = 0) 
as a guide. Specifically, we consider the resonance case 
A c = u! m . Under the condition (A c + cu m ) ^> xo, 
we discard the terms ±ix(t)e ±l< - Ac+Wm ^ in Eq. ([6]) by 
RWA, then by letting u(t) = {SA*8B) + {SB*SA), v(t) = 
i{{5B l 5A) - (8AUB)), and w{t) = (5A l SA) - (SB'SB), 
we can obtain the Bloch equations 

u(t) = 4>{t)v(t), (14a) 

v(t) = -j>(t)u(t)+2 X (t)w(t), (14b) 
w(t) = -2 X (t)v(t). (14c) 

Under the initial condition Ui — 0, Vi — 0, and Wi — 
—n m , we find that, when 



1 



the solution of the Bloch equations is [46 



(15) 



u(t) = —v(t) — ™^ sech[o;(t — to)]j (16a) 
ol 2xo 

w(t) = n m tanh[a(t — to)]. (16b) 
Therefore the average quasi-phonon number evolves as 

{StfSb) = (SB* SB) = 22(1 - tanh[a(t - to)]). (17) 

When a(t — to) ^> 1, we have (5b*5b) w (for example, 
tanh5 = 0.999909), which implies that thermal noise in 
the mirror can be extracted almost completely. Thus, 
Xo = \ V a 2 + /? 2 is a relation to implement efficient cool- 
ing of the mirror in the absence of dissipation. 

It is useful to note that for the constant-pulse coupling 
case [i.e., x(*) is a constant], one can also extract energy 
from the mirror, but the corresponding solution is oscil- 
latory at a Rabi frequency, i.e., the phonon number in 
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FIG. 3: (Color online) (a) Time evolution of the displaced 
phonon number {Sb^Sb} and displaced photon number (Sa^Sa) 
in the chirped-pulse coupling case for two different values for 
the 7c. The solid and short dashed curves are for the case 
of 7c/oj m = 0.0435, while the dashed and dash-dotted curves 
are for the case of "/ c /^m = 0.00435. (b) Plot of the real 
and imaginal parts of (b(t)) vs the scaled time uj m t. Here 
A c = u) m . 



the mirror is a cosine function of time. In this case the 
timing of the constant pulse is crucial in order to locate 
the instant when the phonon number is minimum. How- 
ever, our scheme does not have such a timing control issue 
because the solution (fT7| indicates that after the pulse 
duration (about 2to), the residual quasiphonon number 
of the mirror changes slowly in time by the behavior of 
the tanh[a(i — to)] function. 

As a remark, we indicate that xo should be much 
smaller than the mirror frequency u> m in order to meet 
the condition of RWA. The violation of RWA would mean 
that the parametric interaction of the form SalSb^ + h.c. 
becomes important, which generally leads to heating of 
the system. Therefore by Eq. (fl"5]). the pulse duration 
characterized by aT 1 cannot be arbitrarily short, i.e., 
a < 2xo <C 4w m . Our numerical calculations (without 
RWA) indicate that the parameters used in Fig. [3] are 
quite sufficient for RWA to be valid. 



C. Dissipative case 

In realistic experiments, the interactions with environ- 
ments will inevitably lead to dissipation of the system. 
In addition, the counter-rotating terms (So) 8tf + h.c.) ig- 
nored in RWA may modify the dynamic process. In what 
follows, we numerically study the cooling process in the 
dissipative case beyond RWA. We consider the realisti- 
cally experimental parameters for the system 34] : uj m 
2tt x 73.5 MHz, j m « 2tt x 1.3 kHz, 7 C « 2tt x 3.2 MHz, 
and g w 2tt x 843.1 Hz. Namely, j m /u m rj 1.768 x 1(T 5 , 
7 c /w m w 0.0435, and g/u m « 1.147 x 10~ 5 . With these 
parameters and the chirped pulse given in Fig.[5J we solve 
Eq. ([9]) numerically. In Fig. [3ta), we plot the time evo- 
lution of the mean displaced particle numbers. We see 
that the (Sb^Sb) decreases rapidly from its initial value 
(n m = 1000) to a relatively small number (about 34) 



when the chirped pulse is applied. At the same time, 
(5cv8a) in the cavity increases rapidly from zero to a 
peak value and then decreases to zero gradually through 
the cavity decay channel. We note that the residual fluc- 
tuations of the mirror are limited by the noise of the 
system, mainly of the cavity-field damping. Our numer- 
ical investigations show that when the cavity-field decay 
rate is 7 c /w m « 0.00435, the residual (db^Sb) can fur- 
ther be reduced to 1.08 [the dash line in Fig. Efa)]. In 
addition, numerical calculations indicate the correction 
from the counter-rotating terms is negligible with these 
parameters. 

We remark that during the pulse interaction, the mir- 
ror attains a non-zero coherent amplitude (i.e., {b(t)) ^ 
0) according to Eq. (|13b[) . But such a coherent motion 
should not be confused with the fluctuations we aim to 
reduce in this paper. For the present parameters used in 
Fig.C2a), we plot the time evolution of (b(t)) in Fig.[3jb). 
We see that there is a small amplitude |(6(2to))| ~ 2 near 
the end of the pulse interaction. Actually, when the ther- 
mal fluctuation of the mirror is completely transferred to 
the cavity, the mirror in the displaced representation will 
be in its ground state. Therefore, in the original repre- 
sentation, the mirror is prepared in a coherent state of 
the mechanical motion. 

We also point out that although there are residual cav- 
ity photons after the pulse duration (around 2to), the 
heating due to such photons is found to be negligible 
because of the weak coupling strength g. We can es- 
timate that for a residual cavity photon number n r at 
time 2io, the mirror can be excited to have phonon num- 
ber (gn r /u! m ) 2 , which is of the order of 0.001 with the 
parameters used in Fig. [3] (-f c /uj m — 0.0435). Therefore 
the heating is mainly due to the heat bath of the mirror. 
For example, we find that the heat bath of the mirror 
would increase the phonon number from 34 to 38 when 
the time evolves from u> m t = 80 to 300. 



D. Effects of the j3 parameter 

The frequency modulation characterized by the param- 
eter j3 is a main feature of the chirped coupling. In the 
case of (3 = 0, the coupling corresponds to a t pulse 
because J_ asech[a(t — to)]dt = ir is the pulse area. 
However, such a simple 7r pulse generally does not bring 
an optimal cooling when dissipation and counter-rotating 
terms are included. The parameter /3 therefore provides a 
way to adjust the pulse for a better cooling performance. 
In Fig. U we demonstrate this feature numerically by 
plotting the final mean displaced phonon number of the 
mirror (defined at time t — 2to) as a function of /3. There 
are two situations [Fig.QJa) and0|b)] that we will discuss 
below, but in both figures, it is apparent that non-zero 
values of (3 can better reduce the displaced phonon num- 
ber of the mirror. 

We point out that the final displaced phonon number 
can become less sensitive to the detuning A c when |/3| 
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FIG. 4: (Color online) (a) Plot of the final mean displaced 
phonon number (8b*8b)f vs the phase modulation amplitude 
f3 for various detunings A c /uj m = 1 and 1.02. (b) Plot of 
(5b* Sb}f vs p for various parameters 8 = —0.1, 0, and 0.1. 



is increased. This is shown in Fig. Eta) where we can 
compare the sideband resonance case A c = u) m with a 
slightly off resonance case. We see that although cooling 
with an off-resonance A c is less effective, the dependence 
on A c becomes weaker as |/3| increases. For the parame- 
ters used in Fig.|4ja), the final displaced phonon numbers 
are essentially the same when |/?| > 0.2w m . This shares 
a similar feature in two-level systems as a chirped pulse 
can make efficient population transfer in the presence of 
inhomogeneous broadening. 

In Fig. B|b), we illustrate the effect of j3 on cooling 
when there are uncertainties in controlling the pulse area. 
Such an error, for example, may come from an inaccurate 
value of the coupling strength g. Let us express x(t) as 



x(t) = (1 + S)xosech[a{t - to)], 



(18) 



with S describing the deviation. At = 0, Fig. @|b) 
shows that a modest change of 5 can affect the final dis- 
placed phonon number quite significantly. In fact, we 
notice that the case 6 = 0.1 in the figure actually cor- 
responds to a better cooling. This indicates that xo in 
Eq. (fl"5l) is no longer optimal for cooling because of dis- 
sipative effects. The search for optimal pulse parameters 
relies on numerical work, but Fig. 0|b) suggests that the 
frequency modulation with a suitable range of j3 may ease 
the sensitivity of 5 and hence improve the cooling perfor- 



mance even though the pulse parameters are not exactly 
optimal. 



V. CONCLUSION AND REMARKS 

To conclude, we have proposed a method to cool a mov- 
ing mirror in cavity optomechanics by a chirped pulse. 
Within the linearization framework, we have shown how 
a chirped pulse coupling can be achieved by an exter- 
nal driving field, and numerically we have demonstrated 
that thermal fluctuations in the mirror can be signifi- 
cantly transferred to the cavity after the pulse. In par- 
ticular, the frequency modulation plays a positive role in 
the cooling process especially when there are uncertain- 
ties in controlling the detuning and pulse areas. 

Finally, we remark that it would be difficult to present 
a general comparison of the cooling efficiency between 
our scheme and the resolved sideband cooling. This is 
because Eq. © has no analytic solution and so the resid- 
ual phonon number can only be calculated numerically. 
Nevertheless, we notice that by decreasing the cavity de- 
cay rate, the residual phonon number can be lowered. 
As a specific example, with 7 c /w TO = 0.001 and the same 
other parameters as in Fig. [3J the residual phonon num- 
ber can reach 0.64. Therefore, the system under such 
parameters may effectively be considered as the ground 
state, although this residual phonon number is hi ghe r 
than the resolved sideband cooling limit (■f c /Aoj m ) 2 {28|. 
The main purpose of this paper is to provide an alterna- 
tive method of cooling based on pulsed interaction, which 
is a transient solution rather than a steady-state one. In 
other words, the process can occur in a finite duration of 
time, and this could be a useful feature for manipulating 
quantum states of the mirror. 
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